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ABSTRACT 



Context. To investigate the distribution of rotation periods of asteroids from different regions of the Solar System and distribution of 

diameters of near-Earth asteroids (NEAs). 

Aims. Verify if nonextensive statistics satisfactorily describes the data. 

Methods. Light curve data was taken from Planetary Database System (PDS) with Rel > 2. Taxonomic class and region of the Solar 

System was also considered. Data of NEA were taken from Minor Planet Center 

Results. The rotation periods of asteroids follow a g-Gaussian with q = 2.6 regardless of taxonomy, diameter or region of the Solar 

System of the object. The distribution of rotation periods is influenced by observational bias. The diameters of NEAs are described 

by a g-exponential with q = 1.3. According to this distribution, there are expected to be 994 ± 30 NEAs with diameters greater than 

1km. 
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1. Introduction 

Asteroids and comets are primordial bodies of the Solar System 
(SS). The study of the physical properties of these objects may 
lead to a better understanding of the processes of formation of 
the SS, and, by inference, of the hundreds of exo-Solar systems 
already known. Distribution of rotation periods and diameters of 
asteroids are two parameters that may give pieces of information 
concerning the evolution of the SS. The first attempt to describe 
histograms of rotation periods of asteroids was made by Harris & 
Burns ( 1979). This work and also others that have followed have 
shown that rotation periods of big asteroids (D >30-40 km) fol- 
low a Maxwellian distribution. Harris & Pravec 



i) have ana- 
lyzed a sample with 984 objects and have confirmed that the dis- 
tribution of rotation periods of asteroids with diameters D > 40 
km is Maxwellian, with 99% of confidence, though this hypoth- 
' esis can be rejected at 95% of confidence. They suggest that 
' objects within this diameter range are primordial, or originated 
from collisions of primordial bodies. It is known that for median 
sized (10 < D < 40 km ) and small (D < 10 km) asteroids, the 
distribution of rotation periods is not Maxwellian. The analysis 
of the data suggests the existence of a spin-barrier for the aster- 
oids with diameters between hundreds of meters and 10 km and 
with more than 11 rotations per day (d ') (period of about 2.2 
h). The absence of a substantial quantity of asteroids with period 
less than 2.2 h may be due to the low degree of internal cohesion 
of these objects. The majority of the sample may contain rubble 
pile asteroids (Davis et al. 1 19791 Harrisll996^ that are composed 
by fragments of rocks kept together by self-gravitation. For ob- 
jects below 0.2 km it is observed rotation periods smaller than 
the spin-barrier, suggesting that these objects have a high inter- 



nal cohesion, implying that they may be monolithic bodies. The 
difficulty in the modelling of rotation periods of asteroids as a 
whole may be associated to the combined action of many mech- 
anisms such as collisions (Paolicchi, Burns & Weidenschilling 
I2002I I. gravitational interactions with planets (Scheeres, Marzari 
& Rossi[2004), angular momentum exchange in binary or multi- 
ple asteroid systems (Scheeres (2002)), or torques induced by so- 
lar radiation, known as YORP effect (from Yarkovsky-O'Keefe- 
Radzievskii-Paddack) (Rubincam l2U00l l. Particularly, YORP ef- 
fect is strongly dependent on the shape and size of the object and 
its distance to the Sun. 

Near-Earth Asteroids (NEAs) is a subgroup of SS asteroids 
whose heliocentric orbits lead them close to the Earth's orbit. 
More than 7000 NEAs are known up to 201 1 . The study of these 
objects is relevant once it may bring information regarding the 
birth and dynamic evolution of the SS. Also a special interest in 
these objects is related to the possibility of collision with the 
Earth with obvious catastrophic consequences (Alvarez et al. 
1980). They also may be potential sources of raw material for 
future space projects. 

The evaluation of the number of asteroids per year that may 
reach the Earth as a function of their diameters is essential for 
the determination of the potential risk of a collision. One of the 
first attempts to estimate this flux was done by Shoemaker et al. 
(fT979l) . 

The impact flux may be taken from the accumulated distri- 
bution of diameters of the NEAs. This distribution is indirectly 
obtained by the current asteroid surveys, according to the abso- 
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lute magnitude H. The distribution of absolute magnitude H is 
described by Jedicke, Larsen & Spahr (I2002I I. 



\ogN ^ aH + p. 



(1) 



where A^ is the number of objects, a is the "slope parameter" and 
j6 is a constant. This relation asymptotically models the observed 
distribution of H. Departure from this power law is probably as- 
sociated with the observational bias due to physical and dynam- 
ical properties of the asteroids (orbital elements, size, albedo), 
and instrumental limitations (CCD, detection software, among 
others). So Eq. ([T]) may be used to describe a given population 
of asteroids if a correction of the bias is made in the raw data. 

The diameters of asteroids D may be given in function of 
their absolute magnitudes and their albedos p^ according to 
Bowell et al. ( fT989] l 



D = 1329- 



10-^/^ 



(2) 



The albedo is the rate of superficial reflection and its value is 
essential for the estimation of the diameters of the asteroids. 
The values of the albedos of asteroids varies according to the 
superficial mineralogical composition (taxonomic complex) and 
object's shape. Typical values range from 0.06 ± 0.02 for low 
albedo objects of C taxonomic complex up to 0.46 + 0.06 for 
high albedo objects of V type (Warner, Harris & PraveC i2009 ). 

Combination of equations ([T]i and (|2]l leads to a power law 
behavior: 



A^(> D) = kD-^. 



(3) 



The parameters are estimated by Stuart (2003'), b - 1.95 {a - 
bis, admitting the same albedo for all the sample) k - 1090, 
and D is given in km. According to this expression and taking 
into account the uncertainties of the measures, Stuart & Binzel 
(I2004I I have estimated that there may exist 1090 ± 180 objects 
with diameters equal to or greater than 1 km {H - 17.8). 

2. Nonextensive statistics 

In order to model the accumulated distribution of periods and di- 
ameters of asteroids, we have applied results from Tsallis nonex- 
tensive statistics. This choice comes from observational evi- 
dences that astrophysical systems are somehow related to nonex- 
tensive behavior. It is known that system with long-range in- 
teractions (typically gravitational systems) are not properly de- 
scribed by Boltzmann-Gibbs statistical mechanics (Landsberg 
[1990 ). Along the last two decades it has been continuously de- 
veloped the nonextensive statistical mechanics that is a gener- 
alization of Boltzmann-Gibbs statistical mechanics. Tsallis pro- 
posed in 1988 (Tsallis il988D a generalization of the entropy. 



So^k 



q-\ ' 



(4) 



where /:>,■ is the probability of the i-th microscopical state, W is 
the number of states, fc is a constant (Boltzmann's constant) and 
q is the entropic index. As ^ — > 1 5 cy is reduced to Boltzmann- 
Gibbs entropy S\ - -kY^f Pi^^ Pi- It was soon realized that the 
nonextensive statistical mechanics could be successfully applied 
to self-gravitating systems: Plastino & Plastino (1993) found 
a possible solution to the problem of the existence of a self- 
gravitating system with total mass, total energy and total en- 
tropy simultaneously finite, within a nonextensive framework. 



Many examples of nonextensivity in astrophysical systems may 
be found. We list some instances. Nonextensivity was observed 
in the analysis of magnetic field at distant heliosphere associated 
to the solar wind observed by Voyager 1 and 2 (Burlaga & Vinas 
12005 . Burlaga & Ness 2009 Burlaga & Ness 2010). The distri- 
bution of stellar rotational velocities in the Pleiades open cluster 
was found to be satisfactorily modelled by a (7-Maxwellian dis- 
tribution (Soares et al. 2 006 ). The problem of Jeans gravitational 
instability was considered according to nonextensive kinetic the- 
ory (Lima, Silva & Santos 120021) . Nonextensive statistical me- 
chanics was also used to describe galaxy clustering processes 
(Wuensche et al. 2004) and temperature fluctuation of the cosmic 
background radiation (Bernui, Tsallis & Villela 120061 Bernui, 
Tsallis & Villela 120071) . Fluxes of cosmic rays can be accurately 
described by distributions that emerge from nonextensive statis- 
tical mechanics (Tsallis, Anjos & Borges 2003 ). A list of more 
instances of applications of nonextensive statistical mechanics 
in astrophysical systems may be found in Tsallis (2009). 

It is important to mention that the index q has a physical in- 
terpretation — it expresses the degree of nonextensivity — and 
for some systems it can be determined a priori (based on dynam- 
ical properties). In fact a nonextensive system is characterized by 
a ^-triplet and not just by a single q (additional information can 
be found in Tsallis (I2009I I). Such triplet was already obtained 
in an astrophysical system (Burlaga & Vinhas 120051 Burlaga & 
Ness 2009 Burlaga & Ness 2010). 

Maximization of S q under proper constraints leads to dis- 
tributions that are generalizations of those that appear within 
Boltzmann-Gibbs context. For instance, if it is required that 
the (generalized) energy of the system is constant (Curado & 
Tsallis II99II 1 then the probability distribution that emerges is a 
^-exponential. 



p{x) OC tX^ (-PqX), 



(5) 



Pq is the Lagrange multiplier (not to confound withjS of Eq. ([T]i)- 
The ^-exponential function is defined as (Tsallis IT994I ) 



[l + (l-^)x]^ 



(6) 



The symbol [d\+ means that [a] = a if a > and [a] = if 
fl < 0. The ^-exponential is a generalization of the exponential 
function, that is recovered if ^ — > 1 . If the constraint imposes that 
the (generalized) variance of the distribution is constant, then 
the distribution that maximizes 5^^ is a ^-Gaussian (Tsallis et al. 
1995, Prato & TsaUis 1999 J. 



p{x) OC exp (-/3,x ). 



(7) 



The g'-Gaussian recovers the usual Gaussian at ^ = 1, and par- 
ticular values of the entropic index q turn p{x) into various 
known distributions, as Lorentz distribution, uniform distribu- 
tion, Dirac's delta (see TsalHs et al. 1995 Prato & Tsaflis 1999] 
for details). 

The Lagrange parameter/?^ also has a precise physical mean- 
ing. Within the statistical mechanics context, the Lagrange pa- 
rameter/?^ in Eq. (|5]l is related to the inverse of the temperature 
(j6i = l/{kT) if q - 1), and in Eq. (|7]i it is related to the in- 
verse of the variance (J3] = I /(2cr^) in normal diffusion, cr^ is 
the variance). Generally speaking, the inverse of the Lagrange 
parameters are associated to the finiteness of the first or second 
moments of the distribution. 

Inverse cumulative distributions of the family of the q- 
exponentials are usually, and conveniently, graphically repre- 
sented by means of log-log plots, and Fig. [T] shows a general 
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Fig. 1. Inverse cumulative distribution of a generalized q- 
exponential, N>{x) - Mexp (-/3qx'>'). j - 1 is a ^-exponential, 
7 = 2 is a ^-Gaussian, M is the size of the sample. Values of the 
parameters for this particular instance are g' = 1 .5, fiq = 0.1 and 
M — lO'*. It is indicated the transition point (x'Y given by Eq. 
(O. Dashed curve (red on line) is an ordinary {q - 1) exponen- 
tial withySi - 0.1, y - 1 and M - lO'*. It is evident that the tail 
behavior is entirely different from both curves. 



instance for the function N>(x) - M exp (-/3qX^). y is a gen- 
eral parameter that recovers the ^-exponential (y = 1), the q- 
Gaussian (y = 2), and in a generalized scheme, it may assume 
other values, i.e. it may represent other distributions. The tails of 
^-exponentials are power laws (exp (-/3qx) ~ [{q - l)yS^x]'^*'"''^ 
for jc » 1 and q > 1) and thus the slope of the asymptotical 
power law regime in a log-log plot leads to the determination of 
the parameter q (slope - y/(l - q) for general y not necessarily 
equal to one). For small values of the independent variable x (we 
are assuming x > 0) this graphical representation displays a re- 
gion that appears to be quasi-flat in this log-log plot; of course it 
is not flat, once the function is monotonically decreasing by con- 
struction. We are calling this region as "quasi-flat" to distinguish 
it from the asymptotical power law region. Figure [1] shows the 
intersection of two straight lines that represent the two regimes. 
The intersection is the transition point between the regimes (the 
crossover), and it is given by 



1 



[iq - DA 



(8) 



Figure [T] also displays an ordinary (^ = 1) exponential, for com- 
parison. Exponentials of negative arguments decay vary fast, 
and when represented in log-log plots this feature becomes 
clear, once the exponential asymptotically presents slope — > -oo. 
Coherently Eq. ([8]) gives jc* — > oo for q = I, indicating that there 
is no crossover of an exponential to a power law regime. 

We have found that the distribution of diameters of NEAs 
follows a ^-exponential and the observed rotation periods of all 
asteroids, regardless of their diameters, mineralogical composi- 
tion or region of the SS, are well approximated by a ^-Gaussian. 
We have used two samples from databases of diff'erent years, in 
order to verify the eff'ect of the observational bias. 

3. Observational data 

One important problem in the evaluation of distributions of rota- 
tion periods and diameters of asteroids (and of course the same 
applies for other observables) is that the data are possibly influ- 
enced by observational bias. In order to take into account this 



eff'ect, we have considered samples from databases of two dif- 
ferent years: 2005 and 2010 for rotation periods, and 2001 and 
2010 for diameters of asteroids. 

Two versions of the Lightcurve Derived Data, available at the 
Planetary Database System (PDS), were used: version 7 (V7), 
with 1971 periods, and version 11 (VI 1), with 4310 periods 
(Harris, Warner & Pravec 120051 12010 1. The periods are classi- 
fied according to a quality code of the reliability of the estimated 
period, defined by Harris & Young ( I1983I I. We have used periods 
with Rel > 2 (Rel from reliability) that means they are accurate 
to =» 20% which resulted in 1621 entries for the V7 and 3567 as- 
teroids for the VI 1. Cross-checking the VI 1 sample with a com- 
pilation of taxonomic classifications, also available at PDS, has 
revealed that about 40% (1487) of these asteroids have approx- 
imately known mineralogical composition. The asteroids have 
been separated into three main classes: C , S and X complexes, 
foflowing the SMASS II system of Bus & Binzel (2002), with 
respectively 503, 663 and 321 objects. The diameters of these 
sub-samples have been calculated with Eq. (|2]l with the absolute 
magnitude H available from MPCORB - Minor Planet Center 
Orbit Database (MPC 2010). 

We have also used two versions of the compilation of abso- 
lute magnitudes H, namely that of Oct., 2001, with 1649 NEA's 
(similar to Stuart's (2001 ) procedure) and Oct., 2010, with 7078 
objects. (MPC 2010). We have adopted p^ = 0.14 + 0.02 for the 
NEAs population albedo. This value was estimated by Stuart & 
Binzel ( 120041 1 and it takes into account the great variety of tax- 
onomic types that are found in the NEAs (Binzel et al. 2004] l. 
In order to estimate the validity of the diameters estimated by 
Eq. (|2]), we have considered the diameters of 101 asteroids ob- 
tained from Spitzer Space Telescope data (Trilling et al. 201(31. 
This resulted in about 20% of error. We considered that this 
value, though not small, is reasonable for the purposes of our 
study. 

4. Distribution of rotation periods 

Fig. |2] shows the decreasing cumulative distribution of peri- 
ods of V7 and VI 1, and superposed (7-Gaussians iN>{p) = 
M exp (-/3q p^)) and it is seen that these functions quite sat- 
isfactorily describe almost all the data with q = 2.0 ±0.1, 
/3q = 0.016 ± 0.001 h'- and M = 1621 (M is flie number of 
objects) for V7, and ^ = 2.6 ± 0.2, /?,, = 0.025 + 0.002 h^^ and 
M = 3567 for VI 1. Parameters were found by a nonlinear least 
square method. Fig.|2]also presents two ordinary Gaussians, and 
it is evident that these q = I distributions are completely unable 
to represent the data. Confidence level for both fits is 95%, ac- 
cording to x^ test. This suggests that the distribution does not 
depend on (i) the diameters; (ii) the mineralogical composition; 
(iii) the region of the SS in which the object is found. The latter 
is particularly important once the sample includes NEAs, trans- 
netunian objects (TNO), asteroids from the main belt (MBA), 
Jupiter Trojans (JT) and dwarf planets like Ceres and Pluto. The 
values of the entropic indexes (q = 2.0 for V7 and q = 2.6 
for VI 1) — rather distant from unit, that is, distant from the 
Maxwellian distribution — may indicate that long-range inter- 
actions play an essential role in the distribution of rotation peri- 
ods. According to Eq. ^ (with x* = p*, y - 2), the transition 
point is p' = 7.91 ± 0.01 h (/ ~ 3 d"') for the data of V7, and 
p* = 5.00 ± 0.02 h (/ ~ 5 d-i) for the data of VI 1. The tran- 
sition points for both samples differ from the critical period of 
the spin barrier, and thus the transition is not a consequence of 
physical processes. Warner & Harris ^20101 ) have demonstrated 
that the periods are more accurately determined for objects with 
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Fig. 2. Decreasing cumulative distribution of periods of V7 
(green dots on line) and VI 1 (black dots on line) of PDS 
(NASA) with Rel > 2, and superposed ^-Gaussians {N>(p) - 
Mexp^{-/3gp^)). V7: q = 2.0, yS^ = 0.0161 h-^, M = 1621; VI 1: 
q = 2.6, y6^ = 0.025 h^^ M = 3567. Fittings of the periods of 
V7 for P > 50 h are not good. This does not happen with VI 1, 
and it possibly indicates the increase in the accuracy of the data. 
Dashed (violet on line) and dot-dashed (magenta on line) lines 
are usual (q - 1) Gaussians, withjSi - 0.0161 h"-, M - 1621, 
andySi = 0.025 h^^, M = 3567. 
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Fig. 3. Log-log plot of the decreasing cumulative distribution 
of periods of 3567 asteroids (dots) with Rel > 2 taken from 
the PDS (NASA) and a ^-Gaussian distribution (N>(p) - 
Mexpg(-/3gp^)) (solid line), with q = 2.6, fig = 0.025 h^^ 
M - 3567. The other curves are 663 S-complex asteroids (di- 
amonds, blue on-line), 503 C-complex asteroids (squares, green 
on-line), 321 X-complex asteroids (triangles, magenta on-line). 
Inset shows the 3567 asteroids and the ^-Gaussian in a linear- 
linear plot. 



periods /:> < 8 h and light curve amplitudes A > 0.3 mag, so we 
conclude that the difference between the transition points of the 
two versions is due to the action of observational bias. 

We have taken separately the taxonomic complexes C, S, 
and X (VI 1), shown in Fig. [3] and we have found that all of 
them are properly described by ^-Gaussians within 95% of con- 
fidence level iq = 2.6 + 0.2 andyS, = 0.021 ± 0.002 h^ for S, 
^ = 2.0 ± 0.1 and fig = 0.015 + 0.001 h^^ for C and ^ = 2.0 ± 0.1 
andfig - 0.010 + 0.007 h"^ for X; fitted curves are not indicated 
in Fig. [3]). The diff'erence between the parameters for each com- 
plex is related to the size of the sample, and they are statistically 
similar The fittings for the C and X complexes can be improved 
if we admit that the number of objects are 10% higher than that 
found in the sample. The transition point p* gets smaller from 
V7 to V 1 1 , and this may be an indication that the fraction of fast 
rotators (/ > 5 d"') may still be sub-estimated. 



5. Distribution of diameters of near-Earth asteroids 

We have found that the decreasing cumulative distribution of 
diameters of NBAs can be fitted by ^-exponentials. The fitting 
of a ^-exponential to the diameters of 7078 NBAs {N>(D) - 
M exp J-figD)), shown in Fig. ID is quite good for the entire 
range of the data, with a confidence level of 95%:^= 1.3 + 0.1 
and fig = 3.0 ± 0.2 km"' (found by nonlinear least squares 
method). This distribution, however, is influenced by observa- 
tional bias. The ^-exponential distribution can be used to de- 
termine the point in which the sample is supposed to be com- 
plete. Figure |4] compares g'-exponentials that fit observed dis- 
tribution of diameters of known NBAs in October of 2001 and 
October of 2010. The observed distribution of diameters of the 
2011 database follows a ^-exponential with q = 1.3 ± 0.1 and 
fig - 1.5 + 0.1 km ' and the same confidence level. The Figure 
also shows usual (q = I) exponentials, and it can be promptly 
verified their inadequacy in the representation of the whole range 
of the data. 

Once the value of q is the same for both 2001 and 2010 sam- 
ples, we may argue that this parameter is not influenced by the 




D(km) 



Fig. 4. Decreasing cumulative distribution of diameters of 
known NBAs in 2001 (1649 objects, green dots) and in 2010 
(7078 objects, black dots). Solid lines are best fits of q- 
exponentials (N>(D) = M exp^i-figD)). Blue Une (2001): q = 

1.3, fig = 1.5 km-i, M = 1649, red line (2010): q = 1.3, 
fig - 3) km"', M - 7078. Usual exponentials (^ = 1) are dis- 
played in the main panel for comparison (dashed violet, with 
fix - 1.5 km"', M - 1649, and dot-dashed magenta, with 
fix =3km-',M = 7078). 



bias in this case, and it reflects real physical processes. Both 
curves are practically identical in the power-law region, and the 
point of transition to the quasi-flat region diff'er, as expressed by 
the different values of fig. The value of q - 1.3 (diff'erent from 
one) is an indication that not only collisional processes are im- 
plied in the formation of these objects. Other mechanisms may 
also be present: the YORP effect may lead to the decrease of the 
period of rotation up to the point of rupture. This fragmentation 
process may yield the formation of binary or multiple systems. 
About (15 ± 4)% of NBAs with D > 0.3 km and rotation pe- 
riods between 2 and 3 h possibly are binary systems (Pravec, 
Harris & Warner 2007). The transition points according to Bq. 
(O (with X* = D\y^ 1) are D* = 2.22 ± 0.05 km (2001), and 
D* = 1.11 ±0.05 km (2010). 
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The sample is complete up to the upper limit of these in- 
tervals, 2.22 + 0.05 = 2.27 km {H = 16) for 2001 basis and 
1.11 + 0.05 = 1.16 km (// = 17.5) for 2010 basis. The num- 
ber of NEA with D > 2.27 km is virtually the same for the 
sample of 2001 and 2010 (166 ± 8 objects). This is a confirma- 
tion of the completeness of the sample up to // ~ 15 (Jedicke, 
Larsen & Spahr, 2002), and the extension of this limit up to 
// ~ 16 (Harris, 2008 ). Once there has been an increase in the 
efficiency of detection and in the number of surveys (Stokes, 
Evans & Larson, '2002', Larson '2007"), we conclude that the pa- 
rameter's^ indicates the limit of completeness of the sample. For 
D > 1.16 km, the 2010 data are best described by a power-law. 
We have found for Eq. ^,k ^ 994 + 30 and b = 2.24 ± 0.01, 
with correlation coefficient R- - 0.987. The value of b corre- 
sponds to a = 0.448 + 0.002 in Eq. ([TJ, that is a reasonable value 
if compared to the slope of 0.44 found by Zavodny et al. (2QQ%. 
The value of q may be found from b using q - \ + lib, thus 
q - 1.446 ± 0.001, that is within the interval found in the whole 
sample. According to the distribution we have found 994 + 30 
asteroids with D > 1 km {H < 17.7), in a close agreement with 
Mainzer et al. ( 1201 11 1. Analysis of distributions of diameters of 
MBA and TNO according to lines similar to those used in this 
work are welcome. 
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ABSTRACT 



Context. To investigate the distribution of periods of rotation of asteroids of different regions of the Solar System and distribution of 

diameters of near-Earth asteroids (NEAs). 

Aims. Verify if nonextensive statistics satisfactorily describes the data. 

Methods. Data was taken from Planetary Database System (PDS) with Rel > 2 (total of 3567 asteroids). Taxonomic class and region 

of the Solar System was also considered. Data for NEA was taken from Minor Planet Center (total of 7078 NEAs). 

Results. The periods of rotation of asteroids follow a ^-gaussian with q = 2.6 regardless of taxonomy, diameter or region of the Solar 

System of the object. The diameters of NEAs are described by a ^-exponential with q = 1.3. According to this distribution, there are 

expected to be 1240 ±180 NEAs with diameters greater than 1 km. 

Key words, asteroids - rotation periods - diameters - nonextensivity 



1. Introduction 

, Asteroids and comets are primordial bodies of the Solar System 
■ (SS). The study of the physical properties of these objects may 
' lead to a better understanding of the processes of formation of 
the SS, and, by inference, of the hundreds of exo-Solar systems 
already known. Distribution of rotation periods and diameters 
] of asteroids are two parameters that may give peaces of infor- 
mation concerning the collisional evolution of the SS. Pravec 
& Harris (I2000I I have analyzed a sample with 984 objects and 
have argued that the distribution of rotation periods of asteroids 
with diameters D > 40 km is maxwellian, with 99% of confi- 
' dence. They suggest that objects within this diameter range are 
primordial, or originated from collisions of primordial bodies. It 
is known that for small and very small asteroids (D < 10 km), 
the distribution of rotation periods is not maxwellian. The anal- 
ysis of the data suggests the existence of a spin-barrier for the 
' asteroids with diameters between hundreds of meters and 10 km 
' and with less than 1 1 rotations per day (period of about 2.2 h) 
(Pravec, Harris & Warner 2007). The absence of a substantial 
quantity of asteroids with period less than 2.2 h may be due to 
the low degree of internal cohesion of these objects. The ma- 
jority of the sample may contain rubble pile asteroids (Davis et 
al. 1 19791 Harris 1996) that are composed by fragments of rocks 
kept together by self-gravitation. For objects below 0.2 km it 
is observed rotation periods smaller than the spin-barrier, sug- 
gesting that these objects have a high internal cohesion, imply- 
ing that they may be monolithic bodies. The difficulty in the 
modelling of periods of rotation of asteroids as a whole may 
be associated to the combined action of many mechanisms such 
as collisions (Paolicchi, Burns & Weidenschilling 2002), grav- 
itational interactions with planets (Scheeres, Marzari & Rossi 



I2004 i). angular momentum exchange in binary or multiple aster- 
oid systems (Bottke et al. 2002), or torques induced by solar 
radiation, known as YORP effect (from Yarkovsky-O'Keefe- 
Radzievskii-Paddack) (Rubincam ,2000j . Particularly, YORP ef- 
fect is strongly dependent on the shape and size of the object and 
its distance to the Sun. 

Near-Earth Asteroids (NEAs) is a subgroup of SS asteroids 
whose heliocentric orbits lead them close to the Earth's orbit. 
More than 7000 NEAs are known up to 201 1 . The study of these 
objects is relevant once it may bring information regarding the 
birth and dynamic evolution of the SS. Also a special interest in 
these objects is related to the possibility of collision with the 
Earth with obvious catastrophic consequences (Alvarez et al. 
1980). They also may be potential sources of raw material for 
future space projects. 

The evaluation of the number of asteroids per year that may 
reach the Earth as a function of their diameters is essential for 
the determination of the potential risk of a collision. One of the 
first attempts to estimate this flux was done by Shoemaker et al. 
(,1979,) . 

The impact flux may be taken from the accumulated distri- 
bution of diameters of the NEAs. This distribution is indirectly 
obtained by the current asteroid surveys, according to the abso- 
lute magnitude H. The distribution of absolute magnitude H is 
described by Michelson fl2004D . 



logN = aH+/3, 



(1) 



where A^ is the number of objects and a and /3 are constants. 
The diameters of asteroids D may be given in function of their 
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absolute magnitudes and their albedos p^ according to Bowell et 
al. (fT989l l 



D = 1329- 



10^ 

Vp7 ' 



(2) 



The albedo is the rate of superficial reflection and its value is 
essential for the estimation of the diameters of the asteroids. 
The values of the albedos of asteroids varies according to the 
superficial mineralogical composition (taxonomic complex) and 
object's shape. Typical values range from 0.06 ± 0.02 for low 
albedo objects of C taxonomic complex up to 0.46 ± 0.06 for 
high albedo objects of V type (Warner, HaiTis & Pravec 2009 ). 

Combination of equations ([TJ and (|2]i leads to a power law 
behavior: 



N{> D) = kD 



-b 



(3) 



The parameters are estimated by Stuart (2003"), b - 1.95, k - 
1090, and D is given in km. According to this expression and 
taking into account the uncertainties of the measures, Stuart & 
Binzel (2004) estimated that there may exist 1090 + 180 objects 
with diameters equal to or greater than 1 km. 

2. Nonextensive statistics 

In order to model the accumulated distribution of periods and di- 
ameters of asteroids, we have applied results from Tsallis nonex- 
tensive statistics. This choice comes from observational evi- 
dences that astrophysical systems are somehow related to nonex- 
tensive behavior It is known that system with long-range in- 
teractions (typically gravitational systems) are not properly de- 
scribed by Boltzmann-Gibbs statistical mechanics (Landsberg 
1 19901 1. Along the last two decades it has been continuously de- 
veloped the nonextensive statistical mechanics that is a gener- 
alization of Boltzmann-Gibbs statistical mechanics. Tsallis pro- 
posed in 1988 (TsaIhs lT988l l a generalization of the entropy. 



w „q 



Sa^k 



q-l 



(4) 



where pi is the probability of the i-th microscopical state, W is 
the number of states, A; is a constant (Boltzmann's constant) and 
q is the entropic index. As ^ — > 1 5 ^y is reduced to Boltzmann- 
Gibbs entropy S ] = -kYjf Pi In Pi- It was soon realized that the 
nonextensive statistical mechanics could be successfully applied 
to self-gravitating systems: Plastino & Plastino (1993) found 
a possible solution to the problem of the existence of a self- 
gravitating system with total mass, total energy and total entropy 
simultaneously finite, within a nonextensive framework. Many 
examples of nonextensivity in astrophysical systems may be 
found. We list some instances. Nonextensivity was observed in 
the analysis of magnetic field at distant heliosphere associated to 
the solar wind observed by Voyager 1 and 2 (Burlaga & Vinhas 
120051 Burlaga & Ness l20U9l Burlaga & NesslMll. The distri- 
bution of stellar rotational velocities in the Pleiades open cluster 
was found to be satisfactorily modelled by a ^-maxwellian dis- 
tribution (Soares et aI. I2006i l. The problem of Jeans gravitational 
instability was considered according to nonextensive kinetic the- 
ory (Lima, Silva & Santos 2002). Nonextensive statistical me- 
chanics was also used to describe galaxy clustering processes 
(Wuensche et al. 2004 ) and temperature fluctuation of the cosmic 
background radiation (Bernui, Tsallis & Villela 120061 Bernui, 
Tsallis & Villela l2007l ). Fluxes of cosmic rays can be accurately 



described by distributions that emerge from nonextensive statis- 
tical mechanics (Tsallis, Anjos & Borges '2003). A list of more 
instances of applications of nonextensive statistical mechanics 
in astrophysical systems may be found in Tsallis ( 120091) . 

Maximization of Sg under proper constraints leads to a dis- 
tributions that are generalizations of those that appear within 
Boltzmann-Gibbs context. For instance, if it is required that 
the (generalized) energy of the system is constant (Curado & 
Tsallis I199II ) then the probability distribution that emerges is a 
^-exponential, 

p(x) oc exp^(-/3^x), (5) 

/3q is the Lagrange multipUer. The ^'-exponential function is de- 
fined as (Tsanis [T994l ) 



exp x= [1 -l-(l -q)x]', 



(6) 



The symbol [a]+ means that [a] = a if a > and [a] = if 
fl < 0. The ^'-exponential is a generalization of the exponential 
function, that is recovered if ^ — » 1 . If the constraint imposes that 
the (generalized) variance of the distribution is constant, then the 
distribution that maximizes Sq is the g'-gaussian (Tsallis et al. 
\T995\ Prato & Tsallis (19991 ) . 

p(x) oc expg(-PgX^). (7) 

The (7-gaussian recovers the usual gaussian at q - 1, and par- 
ticular values of the entropic index q turn p{x) into various 
known distributions, as Lorentz distribution, uniform distribu- 
tion, Dirac's delta (see TsaUis et al. 119951 Prato & TsaUis [19991 
for details). 

We have found that the distribution of diameters of NEAs 
follows a ^'-exponential and the observed rotation periods of all 
asteroids, regardless of their diameters, mineralogical composi- 
tion or region of the SS, are well approximated by a ^-gaussian. 

3. Observational data 

In order to obtain a distribution of asteroids rotation periods we 
used the compilation of 4127 periods available at the Planetary 
Database System (PDS) (Harris, Warner & Pravec 120101 . In this 
sample the periods are classified according to a quality code of 
the rehability of the estimated period, defined by Harris & Young 
(119831 ). We have used periods with Rel > 2 (Rel from reliabil- 
ity) that means they are accurate to x 20% which resulted in 
3567 asteroids taken into account. Cross-checking the samples 
with a compilation of taxonomic classifications, also available 
at PDS, has revealed that about 40% of these asteroids have ap- 
proximately known mineralogical composition. The 1487 aster- 
oids with taxonomic classification have been separated into three 
main classes: C , S and X complexes following the SMASS II 
system of Bus & Binzel (120021 ). with respectively 503, 663 and 
321 objects. The diameters of these sub-samples have been cal- 
culated with Eq. ^ with the absolute magnitude H available 
from MPCORB - Minor Planet Center Orbit Database (MPC 
.2010J . 

Absolute magnitudes of 7078 NEAs were obtained from 
MPCORB. We have adopted p,, = 0.14+0.02 for the NEAs pop- 
ulation albedo. This value was found by Stuart & Binzel (2004]) 
and it takes into account the great variety of taxonomic groups 
that are found in the NEAs (Binzel et al. 120041). In order to esti- 
mate the validity of the diameters estimated by Eq. (|2]), we have 
considered the diameters of 101 asteroids obtained from Spitzer 
Space Telescope data (Trilling et al. 2010). This resulted in about 
20% of error. We considered that this value, though not small, is 
reasonable for the purposes of our study. 
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Fig. 1. Log-log plot of the decreasing cumulative distribution of 
periods of 3567 asteroids (dots) with Rel > 2 taken from the PDS 
(NASA) and a (7-gaussian distribution (solid line), Eq. (|7]l, with 
q = 2.6, yS = 0.028 h"^ and the normalizing constant M - 3567. 
Particular asteroids are indicated with blue dots (see Tables [1] 
and IS . 
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Fig. 2. Log-log plot of the cumulative distribution of periods 
of 3567 asteroids (circles) with Rel > 2, and 663 S-complex 
asteroids (diamonds, blue on-line), 503 C-complex asteroids 
(squares, green on-line), 321 X-complex asteroids (triangles, 
magenta on-line). ^-Gaussian of the fitting for all asteroids is 
also displayed (red curve). 



Table 1. Asteroids separated by the region in the Solar System. 
Data from PDS (Harris, Warner & Pravec l20T0b . 



Object 


Period (h) 


Region 


1999 CU3 


3.7829 


NEA 


Hektor 


6.294 


JT 


Ceres 


9.074170 


MB (dwarf planet) 


Sedna 


10.273 


TNO 


Pluto 


153.2935 


TNO (dwarf planet) 



Table 2. Small fast rotators NBAs at the approximately flat re- 
gion of Fig.[T] Data from Binzel et al. (120021 ). 



Object 


H 


Diameter (km) 


Period (h) 


1998 KY26 


25.6 


0.04 


0.1784 


2000 D08 


24.8 


0.04 


0.022 


1999 SFIO 


24.0 


0.06 


0.0411 


1999 TY2 


23.1 


0.08 


0.1213 



4. Distribution of periods 

Fig. [U shows the decreasing cumulative distribution of 3567 pe- 
riods and a superposed ^-gaussian iN>(h) - M exp (-/3qX^)) and 
it is seen that it quite satisfactorily describes almost all the data 
with q = 2.6 and/? = 0.028 h^^ (^ ^ 3557 j^ the number of ob- 
jects). This suggests that the distribution does not depend on (i) 
the diameters; (ii) the mineralogical composition; (iii) the region 
of the SS in which the object is found. The latter is particularly 
important once the sample includes NBAs, trans-netunian ob- 
jects (TNO), asteroids from the main belt (MB), Jupiter Trojans 
(JT) and dwarf planets like Ceres and Pluto. Some asteroids are 
identified in the Figure and in Tables [T] and |2] The value of the 
entropic index q = 2.6 — rather distant from unit, that is, distant 
from the Maxwellian distribution — indicates that long-range 
interactions play an essential role in the distribution of periods 
of rotation. 

The spin-barrier period of 2.2 h is close to the intersection 
of the asymptotic power-law and the approximately flat region 
of Fig. [T] This result suggests that in the approximately flat re- 
gion we may find monohthic fast-rotating asteroids (Whiteley, 
Hergenrother & Tholen l2002l l with absolute magnitude H > 23 
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Fig. 3. Decreasing cumulative distribution of diameters of NBAs 
(dots) in log-log plot. Red curve is a ^ exponential, Bq. (|5]l, with 
g. = 1.3,yg = 3 km^', and M = 7078 is the number of NBAs. 



and D < 0.2 km. Some of these asteroids are identified in Fig.[T] 
and in Table |2] Consequently in the asymptotic power-law we 
probably find gravitational bound objects. 

The distributions for each complex taken separately (C, S, 
X) is shown in Fig. |2] They can be satisfactorily approximated 
by (7-gaussians (q - 2.7 for S, ^ = 2.0 for C and q = 2.0 for X). 
For C and X complexes, the fittings can be improved if we admit 
that the number of objects are 10% higher than that found in the 
sample. 

5. Distribution of Diameters of near-Earth asteroids 

The decreasing cumulative distribution of diameters of NBAs is 
not described by Bq. ^ once it is only asymptotically a power 
law. We have tried to fit a ^-exponential to the diameters of 7078 
NBAs (N>(D) = M expg(-/3gD)) and the result, shown in Fig. [3] 
is quite good for the entire range of the data. 

The value of ^ = 1.3 (diff'erent from one) is an indication 
that not only collisional processes are implied in the formation of 
these objects. Other mechanisms may also be present: the YORP 
effect may lead to the decrease of the period of rotation up to 
the point of rupture. This fragmentation process may yield the 
formation of binary or multiple systems. About (15 ± 4)% of 
NBAs with D > 0.3 km and periods of rotation between 2 and 
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3 h possibly are binary systems (Pravec, Harris & Warner |2007| |. Wuensche, C. A., Ribeiro, A. L. B,. Ramos, F. M., Rosa, R. R. 2004, Physica A, 

According to the distribution we have found there are 1240+ 180 344, 743 

NEAs with D > 1 km, that is superior than the 1090 asteroids 
found by Stuart & Binzel (2004). The analysis of the distribu- 
tions of diameters of MBA and TNO are welcome. 
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